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£SJ . Special stochastic representation of the wave function in Quantum 

> | ■ Mechanics (QM), based on soliton realization of extended particles, is 

' suggested with the aim to model quantum states via classical computer. 

Entangled solitons construction being introduced in the nonlinear spinor 
field model, the Einstein — Podolsky — Rosen (EPR) spin correlation is cal- 
, culated and shown to coincide with the quantum mechanical one for the 

spin-1 /2 particles in the singlet state. The concept of stochastic qubits is 
used for quantum computing modelling. 
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1 Geometric Quantum Mechanics and its 

Oh! Stochastic Representation 

+-> ■ 

Sh ' Recent years a very fascinating idea to put QM into geometric language attracts 

the attention of many physicists [1]. The starting point for such an approach 
qh' is the projective interpretation of the Hilbert space TL as the space of rays. 

To illustrate the main idea it is convenient to decompose the Hermitian inner 
product (-|-) in TL into real and imaginary parts by putting for the two L%- 
vectors | -01 ) = U\ + iv\ and \ip 2 ) = u 2 + iv 2 : 

(ihlih) = Gtyi,^)-*n W-i, *h), (i) 

where G is a Riemannian inner product on TL and ft is a symplectic form, that 
is 

G [tpi,ip2) = (ui,u 2 ) + (vi,v 2 ); fl (ipi,ip 2 ) = (vi,u 2 ) - (ui,v 2 ), 

with (•, •) denoting standard L 2 inner product. The symplectic form revealed 
in (1) can acquire its dynamical content if one uses the special stochastic rep- 
resentation of QM suggested in [2, 3, 4]. 

As a first motivation for introducing stochastic representation of the wave 
function let us consider the de Broglie plane wave 

tp = A e-' lkx = A e -^*+<kr) ( 2 ) 
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for a free particle with the energy u>, momentum k, and mass to, when the 
relativistic relation 

,2 2 i 2 2 

K = U) — K = TO 

holds (in natural units H = c= 1). 

Suppose, following L. de Broglie [5] and A. Einstein [6], that the structure of 
the particle is described by a regular bounded function u(t, r), called hereafter as 
soliton, which is supposed to satisfy some nonlinear equation with the Klein — 
Gordon linear part. Let £q = 1/to be the characteristic size of the soliton 
solution u(t,r) moving with the velocity v = k/w. 

Now it is worth-while to underline the remarkable fact behind this re- 
search [7], namely, the possibility to represent the de Broglie wave ^ as the 
sum of solitons located at nodes of a cubic lattice with the spacing a ^> £q: 

Ae-' fa =^ U (f,r + d), (3) 

d 

where d marks the positions of lattice nodes. To show the validity of (3) one 
can take into account the asymptotic behavior of the soliton in its tail region: 



u{x) = J d 4 ke- tkx g(k)S(k 2 - m 2 ) 
and then use the well-known formula 



implying that 

.4 



2tt\ 3 g(m) 



a J 2m 



The formula (3) gives a simple illustration of the wave — particle dualism, show- 
ing that the de Broglie wave characterizes the assemblage of particles — solitons. 

The shortest way to get the stochastic representation of quantum mechanics 
is modify the formula (3). This can be easily performed if one admits that 
the locations of solitons' centers are not regular nodes of the cubic lattice but 
some randomly chosen points. To realize this prescription, suppose that the real 
field <p describes n particles — solitons and has the form 



k=l 



0(t,r) = »W(t,r), ,4) 

where 



supp n supp 0( fc ') = 0, k^k', 
and the same for the conjugate momenta 



7r(t,r) =d£/d^ = £y fc )(i,r), fa =dcj>/dt, 

where C stands for the Lagrangian density of the field (j>. Let us define the 
auxiliary functions 

^)(t ) r) = -^=K^ fc )+ i 7rW/^) (5) 
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with the constants v k satisfying the normalization condition 



= J d 3 x\^\ 2 . (6) 



Now we define the analog of the wave function in the configurational space 
R 3 " 3 x = {n,...,r„} as 

N n 

* N (t, n, . . . , r„) = (h^N)- 1 ' 2 II ^ fe) (*< r ^)' ( ? ) 

i=i fe=i 

where N » 1 stands for the number of trials (observations) and is the 
one-particle field function for the j-th trial. 
Now we intend to show that the quantity 

PN = T7^ I d 3n x\^ N \\ 



(AV)* 

(AV)"Cl 3 » 

where AV is the elementary volume which is supposed to be much greater 
than the proper volume of the particle ^o 3 = Vo <C AV, plays the role of the 
coordinate probability density. To this end let us calculate the following integral: 



(Avf^E 

(AV) 

where the denotation is used 



f d 3n x \* N \ 2 = (ft" ao- 1 E °« ) . 



tlAv v ' 7 

Taking into account (6) one gets 

(Av) n p N = (^Ny^^AN + S), S = J2 a ij, ( 8 ) 

with AA^ standing for the number of trials for which the centers of particles — 
solitons were located in (Av) n . It is worth-while to remark that due to in- 
dependence of trials and arbitrariness of initial data and, in particular, of the 
phases of the functions <p ^ , one can consider the entities for i ^ j as inde- 
pendent random variables with zero mean values. This fact permits to use the 
Chebyshev's inequality [8] to estimate the probability of the events for which 
\S\ surpasses h n AN: 

P(\S\>h n AN)<(h n ANy 2 (S 2 ). (9) 
On the other hand, in view of trials' independence one gets 

<* 2 >=£<4>- do) 
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Now one can take into account that the wave packets p\ k ^ are effectively over- 
lapped if their centers belong to the proper volume domain Vo- This property 
permits to deduce from (6) and (10) the estimate 

(S 2 ) < a n h 2n -^-V n AN, (11) 
\ / - (Av) n 

where a ~ 1 is the "packing" factor for the nearest neighbors. Inserting (11) 
into (9) one finds the following estimate: 

P {\S\ > h n AN) < (a Vo /Av)™ <C 1. (12) 

Applying the estimate (12) to (8) one can state that with the probability close 
to unity the following relation holds: 

(Av) n p N = AN/N, (13) 

signifying that the construction (7) plays the role of the probability amplitude 
for the coordinate distribution of solitons' centers, with in (13) being the 
corresponding probability density. Now let us consider the measuring procedure 
for some observable A corresponding, due to E. Noether's theorem, to the 
symmetry group generator Ma- For_example, the momentum P is related with 
the generator of space translation Mp = —i\J, the angular momentum L is 
related with the generator of space rotation Ml = J and so on. As a result one 
can represent the classical observable Aj for the j-th trial in the form 

Aj= J d 5 XTTjiM A <j)j=Y j J d 3 Xip* {k) M A k) <pf\ 

k=l 

The corresponding mean value is 

N N n 

j=i ]=ik=i j 

= J d 3n x** N A* N + 0(^), (14) 
where the Hermitian operator A reads 

n 

A = Y,KM A k) - (15) 
fe=i 

Thus, up to the terms of the order Vo/AV <C 1, we obtain the standard quan- 
tum mechanical rule (14) for the calculation of mean values. It is interesting to 
underline that the solitonian scheme in question contains also the well-known 
spin — statistics correlation [4]. Namely, if tp^ is transformed under the rota- 
tion by irreducible representation of 50(3), with the weight J, then the 
transposition of two identical extended particles is equivalent to the relative 
27r-rotation of that gives the multiplication factor (— 1) 2J in ^n. To show 

(k) 

this property, suppose that our particles are identical, i.e. their profiles tpj ' 
may differ in phases only. Therefore, the transposition of the particles with the 
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centers at ri and r 2 means the 7r-rotation of 2-particle configuration around 
the median axis of the central vector line ri — r 2 . However, due to extended 
character of the particles, to restore the initial configuration, one should per- 
form additional proper 7r-rotations of the particles. The latter operation being 
equivalent to the relative 27r-rotation of particles, one concludes that it results in 
aforementioned multiplication of by (— 1) 2J . Under the natural supposition 
that the weight J is related with the spin of particles — solitons, one infers that 
the many-particles wave function (7) should be symmetrical under the trans- 
position of the two identical particles if the spin is integer, but antisymmetrical 
if the spin is half-integer (the Pauli principle). Thus, we conclude that in the 
solitonian scheme the spin — statistics correlation stems from the extended char- 
acter of particles — solitons. However, the particles in quantum mechanics being 
considered as point-like ones, it appears inevitable to include the transpositional 
symmetry of the wave function as the first principle (cf. Hartree — Fock receipt 
for Fermions). Now it is worth-while to discuss the evidence of wave proper- 
ties of particles in solitonian scheme. To verify the fact that solitons can really 
possess wave properties, the gedanken diffraction experiment with individual 
electrons — solitons was realized. Solitons with some velocity were dropped into 
a rectilinear slit, cut in the impermeable screen, and the transverse momentum 
was calculated which they gained while passing the slit, with the width of the 
latter significantly exceeded the size of the soliton. As a result, the picture of 
distribution of the centers of scattered solitons was restored on the registration 
screen, by considering their initial distribution to be uniform over the trans- 
verse coordinate. It was clarified that though the center of each soliton fell 
into a definite place of the registration screen (depending on the initial soliton 
profile and the point of crossing the plane of the slit by the soliton's center), 
the statistical picture in many ways was similar to the well-known diffraction 
distribution in optics, i.e. the Fresnel's picture at short distances from the slit 
and the Fraunhofer's one at large distances [9]. 



2 Random Hilbert space 

As a result we obtain the stochastic realization (7) of the wave function 
which can be considered as an element of the random Hilbert space W ranc j with 
the inner product 

(Vi,V 2 )=M(^V 2 ), (16) 

with M standing for the expectation value. As a rude simplification one can ad- 
mit that the averaging in (16) is taken over random characteristics of particles — 
solitons, such as their positions, velocities, phases, and so on. It is important to 
underline once more that the correspondence with the standard quantum me- 
chanics is retained only in the point-particle limit (AV 3> Vo) for N — > oo. To 
show this [2] one can apply the central limit theorem stating that for N — -> oo the 
wave function \I>jv(i,x) behaves as the Gaussian random field with the variance 

<7 2 =p(t,x), xeR 3 ™, (17) 

where p(t, x) stands for the probability density (partition function) of solitons' 
centers in R 3 ™. Random Hilbert spaces being widely exploited in mathematical 
statistics, for quantum applications they were first used by N. Wiener in [10]. 
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To illustrate the line of Wiener's argument, we recall the general scheme of 
introducing various representations in quantum mechanics. Let \tp) be a state 
vector in the Hilbert space H and A be a self-conjugate operator with the 
spectrum a(A). Then the a-representation is given by the wave function 

ip(a) = (a\tp), 

where 

A\a) = a\a), a e a(A). 

In particular, the famous Schrodinger coordinate g-representation is given by 
the wave function 

v(?) = <<#> = E^>w>, (is) 

n 

with \n) being some complete set of state vectors in H. Wiener considered the 
real Brownian process x(s, a) in the interval [0, 1] 3 s, where a <G [0, 1] is the 
generalized index of the Brownian trajectory and the correlation reads 

l 

J da x(s, a)x(s' , a) = min (s, s'). (19) 
o 

To obtain the quantum mechanical description, Wiener defined the complex 
Brownian process 

z(s\a,0) = -j=[x(s,a)+iy(s,0)]; a,/?e[0,l], (20) 

and using the natural mapping R 3 — > [0, 1], for the particle in M 3 , constructed 
the stochastic representation of the wave function along similar lines as in (18): 

{a,0\1>)= J dz(s\a,p)iP(s), (21) 

s€[0,l] 

with the obvious unitarity property 
l 

ds\lP(s)\ 2 = J J dadi\(n. I\r)\ 2 
o [04] 2 

stemming from (19). 
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3 Entangled solitons and EPR correlations 

In the sequel we shall consider the special case of two-particle configurations 
(n = 2) , corresponding to the singlet state of two spin-1 /2 particles. In quantum 
mechanics these states are described by the spin wave function of the form 

^12 = -7= (I 1 1") ® I 2 i) - |1 I) ® |2 t)) (22) 
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and are known as entangled states. The arrows in (22) signify the projec- 
tions of spin ±1/2 along some fixed direction. In the case of the electrons in 
the famous Stern — Gerlach experiment this direction is determined by that of 
an external magnetic field. If one chooses two different Stern — Gerlach devices, 
with the directions a and b of the magnetic fields, denoted by the unit vectors 
a and b respectively, one can measure the correlation of spins of the two elec- 
trons by projecting the spin of the first electron on a and the second one on 
b. Quantum mechanics gives for the spin correlation function the well-known 
expression 

P(a, b) = ^+((73) ® (trt>)Vi2, (23) 

where a stands for the vector of Pauli matrices cr,, i = 1, 2, 3. Putting (22) into 
(23), one easily gets 

P(a, b) = -(ab). (24) 

The formula (24) characterizes the spin correlation in the Einstein — Podolsky — 
Rosen entangled singlet states and is known as the EPR-correlation. As was 
shown by J. Bell [11], the correlation (24) can be used as an efficient criterium 
for distinguishing the models with the local (point-like) hidden variables from 
those with the nonlocal ones. Namely, for the local-hidden- variables theories the 
EPR-correlation (24) is broken. It would be interesting to check the solitonian 
model shortly described in the beforehand points by applying to it the EPR- 
correlation criterium. To this end let us first describe the spin-1/2 particles as 
solitons in the nonlinear spinor model of Heisenberg — Ivanenko type considered 
in the works [12, 13]. The soliton in question is described by the relativistic 
4-spinor field ip of stationary type 
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u 

V 



satisfying the equation 



(i 7 k d k - C + A(w)) ip = 0, 



(25) 



(26) 



where u and v denote 2-spinors, k runs Minkowsky space indices 0, 1, 2, 3; £o 
stands for some characteristic length (the size of the particle — soliton), A is the 
self-coupling constant, Tp = <^ + 7°, 7 fc are the Dirac matrices. The stationary 
solution to the equation (26) can be obtained by separating variables in spherical 
coordinates r, a via the substitution 



1 



u = 



f(r) 



v = 



g{r)a r 



(27) 



where a T — (ar)/r. Inserting (27) into (26) one finds 

A 



— u + tiaxjjv 

c 



^(/ 2 - 5 > = o, 

- c v + i(av)u - e^v + ^ {f - g 2 ) v = 0. 



In view of (27) one gets 

l{(7\/)v = - 







' 1 " 


('- 
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i{a\/)u 



Finally, one derives the following ordinary differential equations for the radial 
functions f(r) and g(r): 

(' + ?«)" G ' + 5 C -«">/. 

-/'=(^+V) 9 +^(/ 2 -9 2 ) 9 . 
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As was shown in the papers [12, 13], these equations admit regular solutions if 
the frequency parameter uo belongs to the interval 

< w < c/£ . (28) 

The behavior of the functions f(r) and g(r) at r — > is as follows: 

g(r)=C ir , f = C 2 , /'-><), 

where C\, Ci denote some integration constants. The behavior of solutions far 
from the center of the soliton, i.e. at r — > oo, is given by the relations: 



A 

f = -e~ vr , 
r 



r 

g = -g, 



where 



v= (^ 2 -c 2 /c 2 ) 1/2 , B = l^+u/c. 



If one chooses the free parameters £o and A of the model to satisfy the normal- 
ization condition (similar to (6)) 



oc 

J d 3 xtp+(p = j drr 2 (f+g 2 ) = h 
o 

then the spin of the soliton reads 

S = y d 3 x tp + Jtp = ^e z , 



(29) 



(30) 



where e z denotes the unit vector along the Z-direction, J stands for the angular 
momentum operator 

J = — i[rv] + ^ c ® cr , (31) 

and do is the unit 2 x 2-matrix. Now let us construct the two-particles singlet 
configuration on the base of the soliton solution (25). First of all, in analogy 
with (22) one constructs the entangled solitons configuration endowed with the 
zero spin: 

(32) 



<Pl2 = 



V2 



where <p\ corresponds to (27) with r = ri, and Lp\ emerges from the above 
solution by the substitution 







' 1 " 




' " 


ri - 









1 
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that corresponds to the opposite projection of spin on the Z-axis. In virtue of 
the orthogonality relation for the states with the opposite spin projections one 
easily derives the following normalization condition for the entangled solitons 
configuration (32): 

J d 3 Xl J d 3 x 2 ft 2 (fii2 = h 2 . (33) 

Now it is not difficult to find the expression for the stochastic wave function (7) 
for the singlet two-solitons state: 

N 

where tp$ corresponds to the entangled soliton configuration in the j— th trial. 
Our next step is the calculation of the spin correlation (23) for the singlet two- 
soliton state. In the light of the fact that the operator a in (23) corresponds to 
the twice angular momentum operator (31) one should calculate the following 
expression: 

P'(a,b)=M J d 3 Xl j d 3 x 2 ->S>+2(J 1 a)® 2(3 2 b)^ N , (35) 

where M stands for the averaging over the random phases of the solitons. In- 
serting (34) and (31) into (35), using the independence of trials j ^ j', and 
taking into account the relations: 

J+(^ T =0, J 3 (y5 T = -<£ T , J-^=lf l , 

J-f l = 0, J3<P l = ~<P l , J+<p l =<p\ 

where J± = J\ ± % J 2 , one easily finds that 

foe 
J drr 2 ( f+g 2 ) 




(36) 



Comparing the correlations (36) and (24) one remarks their coincidence, that is 
the solitonian model satisfies the EPR-correlation critcrium. 



4 Conclusion. Simulation of qubits by proba- 
bilistic bits 

Now we intend to explain how stochastic qubits introduced previously could be 
simulateded by standard probabilistic bits [14]. To this end one should define 
the random phase 3>j for the j-th trial in our system of n solitons — particles. 
Let (r) denote the standard (etalon) profile for the fc-th soliton. The most 
probable position dj k \t) of the fc-th soliton's center for j-th trial can be found 
from the following variational problem: 

J d 3 x V * {k) (t,r) l p^ (r-df } ) -> max, 



9 



thus giving the random phase structure: 

n „ 

^ = ^ arg / d 3 x^ k \t,r)^ (r-df) . (37) 
k=i J 

The random phase (37) can be used for simulating quantum computing via 
generating the following K random dichotomic functions: 

f.(9.) = sign [cos($, +6,)], s = T~K, (38) 

with S being arbitrary fixed phases. Now recall that the qubit is identified with 
the state vector 

M=a|0>+/3|1>, 

corresponding to the superposition of two orthogonal states |0) and |1), as for 
instance, two polarizations of the photon, or two possible spin-1/2 states. It 
is worth-while to compare the standard EPR-correlation (24) with the random 
phases one for the case of n = 2 particles: 

E(/i/ 2 ) = l--|A0|, (39) 

where Ad = 9\ — 62- The similarity of these two functions (39) and (24) of the 
angular variable seems to be a good motivation for the K qubits simulation by 
the dichotomic random functions (38) popularized in the paper [15]. This very 
simple model of stochastic qubits simulation can be employed for simulating Bi- 
photons, EPR states and other entanglement states. We hope that this model 
will be useful for the Shor's and Grover's Quantum Algorithms realization. 



4.1 Initialization of probabilistic bits 

Let us consider N stochastic qubits constructed according to the above algo- 
rithm. Let us call qubit No. 1 the control or signal one. If the signal qubit 
is green then the value assigned will be |0), if it is red the output will be |1). 
Hence, at each output we shall get N initialized integral qubits. 



4.2 Hadamard transform for probabilistic bits 

Let us consider the initial bit 

k) = ^(|o) + |i», 

with equal probabilities P — 1/2 being assigned to the states |0) and |1). Con- 
sider now the Hadamard transform H sending the states to the basis rotated 
over 7r/4: 

ff|o>^-^(|o> + |i», ff|i)^-L(|o>-|i». 

Hence, having applied the Hadamard transform to the bit \q) one gets the final 
bit 

H\q) = \0). 
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4.3 Logical CNOT for probabilistic bits 



The result of application of the logical component CNOT depends on the target 
bit, two possible cases being considered: 

1. If the control bit has the value |1) then the target bit is sent to the opposite 
value. 

2. If the control bit has the value |0) then the value of the target bit is not 
changed. 

All these elementary operations can be realized via classical computer through 
simulating the phase structure of realistic solitons by the generator of random 
numbers connected to the model solitons' generator, e. g. Kerr dielectric with 
the optical excitations or magnetic with the excitations of localized spin inver- 
sion domains. 
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